Transport through single-level systems: spin-dynamics in the nonadiabatic regime 
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We investigate the Fano- Anderson model coupled to a large ensemble of spins under the influence 
of an external magnetic field. The interaction between the two spin systems is treated within a 
meanfield- approach and we assume an anisotropic coupling between these two systems. By using 
a nonadiabatic approach we make no further approximations in the theoretical description of our 
system, apart from the semiclassical treatment. Therewith we can include the short-time dynamics 
as well as the broadening of the energy levels arising due to the coupling to the external electronic 
reservoirs. We study the spin dynamics in the regime of low and high bias. For the infinite bias 
case we compare our results to ones obtained from a simpler rate equation approach, where higher 
order transitions are neglected. We show, that these higher order terms are important in the range 
of low magnetic field. Additionally, we analyze extensively the finite bias regime with methods from 
nonlinear dynamics and discuss the possibility of switching of the large spin. 
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I. INTRODUCTION 

A quantum dot typically consists of 10 5 atoms, elec- 
trons tunneling through such devices experience hy- 
perfine interaction with the nuclear spins of the host 
material^ ^— . Combinations of huge numbers of spins can 
be described as a large external effective spin system 
which interacts with the single electron spin. In exper- 
iments with quantum dots, features like a large Over- 
hauser field 3 - and self-sustained current oscillations have 
been observed^. 

The effects of spin-orbit and hyperfine coupling is in 
the scope of recent research for various kinds of systems, 
like self-assembled quantum dots^£, carbon nanotubes 7 - 
or molecular magnets^—. The latter are promising can- 
didates for spintronicsi^. 

The interaction of electron spins with a large spin 
shows interesting nonlinear effects and the system is 
known to exhibit chaotic behavior— - —. These effects 
appear for the closed system with anisotropic coupling 
and an external magnetic field. Lopez-Moniz and co- 
workers^ discussed the coupling to two external leads, 
which were assumed to be polarized. Within a rate equa- 
tion approach they found, that the chaotic behavior sur- 
vives for small magnetic fields. 

Using a nonadiabatic approach in this paper, we ex- 
tend this approach to the finite bias regime, which is 
not accessible within the rate equation method. Further- 
more, the rate equations method is also restricted to first 
order transitions, which we also extend with our nonadi- 
abatic approach. 

The nonadiabatic approach works well for systems 
whose quantum fluctuations are assumed to be small and 
that are coupled to a nonequilibrium environment. Here, 
we adopt this method (that we had extensively tested for 
nanoelectromechanical systems 20 ) to the new and more 
complex situation of a collective spin instead of an oscil- 
lator variable. This semiclassical description is suitable 
for this systems as long the number of spins, which build 



up this ensemble, is comparatively large2i~— . For a fur- 
ther analysis of correlations between the spin states, for 
instance to study entanglement of the ensemble and the 
single electron spin, a quantum description is certainly re- 
quired. However, whiting a semiclassical description it is 
possible to discuss the main spin dynamics. Even though 
the nuclear spin dynamics can be assumed as slower than 
the electron spin-dynamics 1,13 , we found that in certain 
parameter regimes an adiabatic approximation is not suf- 
ficient to describe the system's dynamics. 

Our paper starts with a detailed description of the 
model and the adaptation of the nonadiabatic approach. 
We directly derive the rate equation approach from our 
nonadiabatic method, which is presented in Sec. Ill B[ Af- 
terwards, we compare the results of both methods in the 
infinite bias regime. In Sec. lIIII we start the investigation 
of the finite bias regime with a dynamical analysis based 
on an adiabatic approach with Green's functions. There- 
with we interpret the nonadiabatic results for the spin 
dynamics. Within our conclusion in Sec HVl we discuss 
the advantages and disadvantages of the methods used 
here. 



II. MODEL 

We assume that a vertical magnetic field B z is applied 
to a Fano Anderson mode l 26 ! 27 . The magnetic field leads 
to a Zeemann splitting of the electronic level Ed into two 
levels 2 ^, corresponding to e a = Ed ± \B Z1 with a Gt, 
Without further interactions this solely leads to two spin- 
dependent current channels. Only electrons with spin- 
up (-down) can tunnel through the upper (lower) energy 
level. These energy levels are broadened due to the cou- 
pling to the leads and for small magnetic fields an overlap 
of both channels exists, but there is no communication 
arranged between the two energy levels and spin-flips 
cannot occur. Here we enable transitions between the 
two energy levels with the help of a large external spin, 
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FIG. 1: Sketch of a single- level system coupled to a large 
external spin J. Here S denotes the electronic spin and B 
the external magnetic field. The detuning of the leads chem- 
ical potentials leads to a transport window in the size of the 
applied bias Vbias- 



which interacts with the electronic spin. Fig. [I] depicts 
a sketch of the considered model. There, S denotes the 
electronic spin operator for the levels, which components 
are defined via 



S. 



s z J- (44-44), 



(i) 



introducing the creation/annihilation operators 4/4 of 
the electronic levels. The vertical magnetic field B z cou- 
ples to the S z operator, leading to the splitting of the ini- 
tial single level. The Hamiltonian for the Fano- Anderson 
model reads 

n FA =^2e d 44 + £ kaacl a(J c kaa 

<7 kaa 

+ £ (^4aX + VLJU ka „) + BJ Z - (2) 

Note, that the left (a = L) and right (a = R) lead oper- 
ators c\, aa /c ka(7 are spin-dependent. This enables us to 
consider polarized leads, where the density of states for 
spin- up and spin-down electrons with energy Skaa differ. 
This could be realized with ferromagnetic leads&22r— , 
leading to a spin-dependent current through the system. 
The third term in the Hamiltonian describes the transi- 
tions between a state in the lead and the electronic levels 
with tunneling amplitude Vkaa- For simplicity, we in- 
clude the prefactor of the last term in Eq. (j2j), containing 
the electronic ^-factor, into the definition of the magnetic 
field. 

As well as the electronic spin operators, the large exter- 
nal spins z-component J z couples to the magnetic field. 
The free motion for it is described by 



where we assume the same g-factor and therewith the 
same magnetic field as for the electronic spin operators. 
This is a simplification in theoretical approaches, cf . 13 i 19 , 
but a generalization is straightforward. 

Here, a large external spin means an effective spin de- 
scribing a big ensemble of spins, for example, the col- 
lective spin of the nuclei in a quantum dot. Electrons 
which tunnel through this device, experience an interac- 
tion with the effective spin of the whole system. This 
interaction is described by Fermi contact hyperfine cou- 
pling termi^ 



V — y ^ Xi SiJi, 



(4) 



B Z J Z 



(3) 



introducing the coupling constant A^. If these coupling 
constants are equal for all components A^ = Aj, one 
speaks of an isotropic coupling. There, the two spins 
decouple in the long time limit. In the anisotropic case, 
where at least for two components A^ ^ Xj is valid, in- 
teresting dynamical behavior was observed^. 

We treat the interaction of the two spins in a semi- 
classical manner. Therefore, we have to assume that 
the quantum fluctuations in the system are small. This 
should be valid as long as the external spin is large and 
its fluctuation can be neglected. As a consequence of this 
assumption, we have no spin decay due to dissipation and 
the large spin's length j is conserved. 

Using a mean-field approximation^ for Eq. (j4]) we ob- 
tain 

v M f = ^2 Xi + J* ^ ~ ^ i = x ^y^ z - 

i 

(5) 

Thereby, the fluctuations SAi = Ai — (Ai),A G S,J have 
been neglected. Now we can build up a closed system of 
equations for the considered system. 

For the large spin we use the commutation relations to 
derive the Heisenberg equations of motion 

(I) = - (A,(4> + B z ) (Jy) + \y(Sy)(J Z ), 

(j y ) = (Xz(S z ) + B z ) (J x ) - \ X (S X )(J Z ), 

(l)= \ x (S X ){Jy)-Xy{Sy){J x ). (6) 

This is a strongly nonlinear system, since the values of 
the electronic spin components depend on the large spin. 

A decisive role in the theoretical description of this sys- 
tem plays the treatment of the electronic spin. By using 
a rate equation approach the electronic spin components 
are obtained from equations of motions similar to Eq. ((6]). 
Therewith the short-time dynamics is included, but the 
contributions of the leads come only in as rates. Follow- 
ing from that higher order transitions are neglected and 
one is restricted to the infinite bias regime. One way of 
including higher order transitions is by applying an adi- 
abatic approximation, where the movement of the large 
spin is assumed to be slow compared to changes in the 
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electronic subsystem. The electronic spin operators can 
than be derived via Keldysh Green's functions^ — and 
are given in explicit expressions. The disadvantage of 
this method is that the short-time dynamics is missing. 

In the next part of this paper we derive equations of 
motion for the electronic spin operators using a com- 
pletely nonadiabatic approach. The latter enables us to 
include the short-time dynamics as well as higher order 
transitions. Following from that we can probe the rate 
equation and the adiabatic approach. 



A. Nonadiabatic approach 

In the framework of the nonadiabatic approach, we cal- 
culate all system quantities by considering the full time- 
dependence of all system quantities. In the following, we 
assume an anisotropic coupling \ y = and A^ = \ z = A, 
and together with the mean-field approach we obtain an 
effective Hamiltonian for the electronic levels 



n c (t) =J2^(t) dUo 



(j x (t)) (4^ 



d ^ d^. 



r{t) = e d ±^-±±(j z {t)), 



a =t,|- 



(7) 



The term Ht contains the coupling to the leads. Based 
on this effective Hamiltonian we can describe the effects 
arising due to the coupling to a large external spin. The 
(J z (t)) - component solely leads to an additional shift of 
the electronic levels 36 , but the coupling to the (J x (t)} - 
component enables transitions between both levels. This 
Hamiltonian corresponds to a two-level or a parallel dou- 
ble dot system. 

The derivation of the equations of motion for the 
spin operators is performed by starting from the Heisen- 
berg equations of motion. The derivation is similar to 
those used before for the description of nanoelectrome- 
chanical systems, for details see Refj2S. We use a flat 
band approximation leading to the spin dependent tun- 
neling rates T aa = 27r^ fc \ V ka a\ 2 5(u ~ Finally, 
the result for the spin operator expectation values yield 

(r = r ff = £ a r a(T ) 
j t (s x (t)) = - r(s x (t)) - (b z + \(J z (t))) (Sy(t)) 



+ ^2 / dujRe [ B n( u > l ) + B t^, t)] , 
(s z + A(J z (i)>) (S x (t)) 



d 
~dt 



T(S y (t)) 
- \(J x (t))(Sz(t)) 

+ W dojlm [B^(lo, t) - Bft(w, t)} , 

a J 

(S z (t)) = -T(S z (t)) + X(J x (t))(S y (t)) 

+ I duRe [B^(u,t) - I%(w,t)] , (8) 



with the definition 

k 

(9) 

for the lead-transition functions 

jB« a {u,t) = ~ i(e a (t)-w- \v)B« a {u,t) 

jBZ a , (w,t)=- i{e a , (t)-u- l -T)B: a , ( W , t) 
.A 



i-(J x {t))B% a {u,t). 



(10) 



The time-dependent z-component of the large spin is in- 
cluded into the effective levels Sait), cf. Eq. ([7j). Here, 
the Fermi function f a {uj) does not depend on the spin 
due to the assumption, that the chemical potentials for 
spin-up and -down electrons in each lead are equal. 

For the numerical calculations we decompose the B^ a 
into their real and imaginary part and discretize the inte- 
gration over the lead energies in N intervals. Therefore, 
we obtain [6 + 16 (TV + 1)] coupled equations, including 
the equations of motion for the large spin, cf. Eq. (|6j). 

Using the definitions for the lead-transitions function 
Eq. (fT0]h the electronic current is obtained from 



J ckv 2 Re [B5 a (w,t)] - T aa (h a (t)) 



(11) 



where e equals the electron charge and the equation of 
motion for the occupation yields 

!<Mt)> =-T(n a (t))±X(J x (t))(S y (t)) 

+ J2 [dio2Re[B^(u,t)}. (12) 

a ^ 

Here, the upper/lower lower sign refers to t / | - elec- 
trons. Note, that the total occupation number of the dot 
(N(t)) = (n^(t)) + (%(£)) still couples to the spin oper- 
ators via the transition functions. In the next section we 
derive a rate equation approach where this dependence 
is omitted. 



B. Infinite bias: rate equation approach 

We recover previous results for this model^ by apply- 
ing an adiabatic approximation to our approach. In con- 
crete terms, we assume the large spin's movement as slow 
compared to the electrons which are entering the system. 
As a consequence, we neglect the time-dependence of the 
large spin in the equations for the lead-transition func- 
tions Eq. ([TO]) , leading to e a (t) = e a and a decoupling 



from the equations for the large spin, (J x (i)) = (J x ). 
The remaining two coupled equation can be solved via a 
Laplace transformation. Note, that this method is not a 
complete adiabatic approach, because in the end we still 
solve equations of motions for the electronic spin compo- 
nents, which is not the case for a full adiabatic ansatz. 
Starting from Eq. ([TO]) we obtain for t —> oo 



-i{e a > 



N(u) 



(13) 



^)-^(4) 2 . Be- 



with N(uo) = (s a > -oj-iT)(£ a -Gj- 

fore inserting these results into the equations for the spin 
operators Eq. J5J), we separate them into real and imag- 
inary part and perform the integrations over uo. Hence, 
the spin operator equations become 

r<&(*)> 



dt 
d 

dt 



dt 



(s x (t)) 

(Sy(t)) 



B z 

(b z 



T{Sy(t)) 

-\{j x (t))(s z (t)), 

(S z (t)) = -T(S z (t)) + X(J x (t))(S y (t)) 

+ 1 (r Lt - r u ) . 



HUt))) (s v (t)), 

A(Jz(t)>) (S x (t)) 



(14) 



This result coincides with 19 : the whole system has now 
been reduced to six coupled equations and the current 
simplifies to 

Jaa(t) = eT aa {S aL - (h a (t))} , (15) 
with the occupation 

f t (na(t)) = - r<M*)> ± X(J x (t))(S y (t)) + T La . (16) 

As mentioned before, the total occupation number of the 
dot (N(t)) decouples from the remaining equations and 
becomes constant in the long-time limit 

(r Lt + r u ) 



(N(t -> oo)> 



(17) 



The equations for the spin operators Eq. and the 
large external spin Eq. (J6j) provide further possibilities for 
an analytic investigation. The fixed points of the system 
can easily be calculated, seei^. 

By varying the magnetic field B z and the tunneling 
rate T three regions of different dynamical behavior were 
obtained from the rate equation approach. In Fig.[2]these 
regions are depicted. The solid black lines correspond 
to parameters used ir*i£. There, the length of the large 
spin equals j = \J\ = 10 and the tunneling rates match 
2IYt = = T. This choice of tunneling rates implies 
Tr| = 0, because in this case the current flows exclusively 
through the upper electronic level. Following from that, 
spin-down electrons get trapped into the lower level and 
only contribute to transport after a spin-flip. 
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FIG. 2: Regions of different dynamical behavior concerning 
the parameters F and B z using the rate equation approach. 
As a complement, the gray lines in Fig. [2] illustrate the change 
of the first two regions if the length of the large spin is varied, 
while the ratio of the tunneling rates are fixed. For a smaller 
spin (j = 10 2), region II decreases to the benefit of region I, 
which increases. Equally, the reverse situation can be seen, if 
the length of the spin is increased (j = 10 2). Modifying the 
ratio of the tunneling rates shifts the vertical line separating 
region II from region III. 



C. Infinite bias: results 

In Fig. [3] the results for regime I are depicted. The 
rate equation results describe damped oscillations as ex- 
pected. For large times the z-component of the large 
spin becomes polarized parallel to the magnetic field and 
one spin-down electron gets trapped in the lower energy 
level. There the spin trajectories end up in one of the 
fixed points 



<j?, 2 > = (o,o,±j), <sr >2 >= 0,0, 



-14 



2T 



(18) 



These points exist in the whole parameter regime and are 
independent of the magnetic field. 

In the nonadiabatic case the dynamical behavior is 
quite different. For small times the spin trajectories fol- 
low the damped results from the rate equation approach. 
But the damping decreases strongly after the time step 
tX ~ 1000. The same amount of time steps later, the am- 
plitude drops down for one oscillation period, followed by 
a return to the initial oscillation. This behavior is similar 
for all spin components. If we consider the z-component 
for the large spin, the turning point appears when it ap- 
proaches its fixed point value (J* x ) = j. 

By varying the parameters we found no damped os- 
cillations in region I at all. For a certain range for the 
magnetic field and the tunneling rate the dynamics ap- 
pear as in Fig.(3j For T/X « 8 the oscillations disappear 
and we enter region II. 

In Fig.2] the results for regime II are depicted. The 
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FIG. 3: Comparison of the nonadiabatic approach and the 
rate equation approach results for regime I. The magnetic 
field equals B z /X = 0.1 and the tunneling rate is chosen as 
Y/X = 9. Initial conditions are (J x (0)) = (J y (0)) = 5(^5 - 
l)/(2>/2), (Jz(0)> = 5/>/2>/5+V5 (S y (0)) = (S z (0)) = 0, 
(S x (0)) = 0.5. 



rate equation approach predicts three kinds of dynami- 
cal behavior, namely strong damped, self-sustained and 
chaotic oscillations. 

For B z /X = 0.2 and small tunneling rate T/A = 0.7 
the system performs strong damped oscillations and runs 
into one of the fixed points 



3-6/ 



0,Si, 2 , 



m 



A 2 



A 



' A 



(19) 



where ((S, J)s4)/((S, J)% q) correspond to 61/62, which 
yield 



14 — J- Lf 



2T 



El 
A 



(20) 



This result coincides perfectly with the nonadiabatic 
result, see left row in Fig.|U The spin components run 
into the fixed point (J3), cf. Eq. ([T9]) . The large spin 
becomes almost completely polarized perpendicular to 
the magnetic field. 

In the first instance, by further increasing the tun- 
neling rate, the rate equations forecast self-sustained 
oscillations for the systems trajectories, followed by a 
chaotic oscillating behavior. Using the nonadiabatic ap- 
proach, things change again. There, the chaotic behav- 
ior appears already in the parameter region, where self- 
sustained oscillations were predicted by the rate equation 
approach. Within the nonadiabatic approach we do not 
recover these self-sustained oscillations in region II, we 
solely observe strong damping or chaotic-like behavior. 
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FIG. 4: Behavior of the spin components in regime II. Left 
graph depicts results for B z /X = 0.2 and T/A = 0.7, where 
fast damping of the trajectories appear in both approaches. 
The middle graph shows results for B z /X = 0.1 and Y/X = 
0.16 where the rate equations predict self-sustained oscilla- 
tions, but for the nonadiabatic results already chaos appear. 
Just as in the right graph, were the tunneling rate is further 
increased (Y/X = 0.015). For clarity, we omitted the rate 
equation results for the chaotic-like regions. The lowest row 
depicts the corresponding current results, for the chaotic-like 
regions the green line denotes — Xr^, the light blue line Z14 
and the dark blue line Z14. Initial conditions are chosen as 
for Fig.03 



For higher values of the external magnetic field (regime 
III) the system performs self-sustained oscillations. The 
results are depicted in Fig. [5] There, the oscillation fre- 
quency of the large spin is close to the Larmor frequency 
uj l — B z /\. In the case without coupling between the 
two spin systems (A = 0), the large spin oscillates ex- 
actly with the frequency ujl in the xy— plane. The z- 
component is fixed due (J z ) = and possesses no cou- 
pling to the magnetic field. 

The left graphs in Fig.[5j depict results for Y = 10/ A. 
Here, the frequency of the electronic spin ^/-component 
matches the frequency of the large spin x- and y- 
component ujs v = = uoj y w 0.8/A. But the z- 

component of the electronic spin is almost twice the fre- 
quency of the other components uos z ~ 1.6 A. This be- 
havior survives for a smaller tunneling rate, right graphs 
in Fig.[5j but there the frequency is even closer to the 
Larmor frequency, uoj x ,j y ,s y ~ 0.96A and uos z ~ 1.91A. 
The appearance of the frequency doubling for (S z (t)) is 
explained in Sec. El 

The electronic current in regime III performs periodic 
oscillations. The results for the rate equation approach 
coincides with the nonadiabatic results, which are de- 
picted in Fig. [6] for the same parameters as in Fig. [5] The 




FIG. 5: Results for regime III, where self-sustained oscil- 
lations appear. In the left graphs (T/A = 10) the oscilla- 
tions are smoothed sinusoidal. For lower tunneling rate, the 
(Si) trajectories perform non-sinusoidal, but periodic oscil- 
lations, depicted in the upper graph on the right side. Pa- 
rameters are B z /X = 1.0, ed/X = and with the initial 
conditions (J x (0)) = (J y (0)) = (J z (0)) = -1 and 

(5 a (0)) = <&(0)> = 0, (S y (0)) = 0.5. 



frequency of the current oscillations uo = 2tt/T matches 
uus z . Not surprising, because this spin component cor- 
responds to the occupation difference of the dot system 
2(S Z ) = (n^) — (ft]). Inserting the latter into the current 
equation Eq. (fT5|) and using the solution for the total oc- 
cupation number in the long-time limit, cf. Eq. (fT7]h we 
obtain for the current through lead a in the rate equation 
frame 



J a (t) =e 



-(s z (tWa t -r a i) 



(21) 



And with T = + Tr^ it is obvious that 2l(£) = 
— 2r(£) for all times t and hence current conservation is 
ensured. The latter is also valid for the nonadiabatic 
approach. 

From the above comparison we can already conclude, 
that higher order terms in the system-lead coupling as 
included in the nonadiabatic approach, do matter in the 
regime of low magnetic field. 



III. DYNAMICS IN THE FINITE BIAS REGIME 

In the last section our investigation was focused on 
the infinite bias regime. There, the interaction between 
the large spin and the electronic system leads to inter- 
esting nonlinear effects. The rate equation approach is 
restricted to this regime of high external bias. With our 




FIG. 6: Nonadiabatic current results corresponding to Fig. [5] 
The dashed-dotted line represents the current for a small tun- 
neling rate, where smoothed oscillations appear. For higher 
T, solid gray line, the oscillation is still periodic but not longer 
sinusoidal. 



nonadiabatic approach we learned, that differences ap- 
pear if one includes higher order transition terms. In 
this section, we take the next step by studying the finite 
bias regime. 

Due to the lack of an easy access to further analytic 
studies for the nonadiabatic approach, we use an adia- 
batic approach based on Keldysh Green's functions for 
the dynamical analysis, for details see Sec.[Bl This en- 
ables us to clarify the effects of a finite transport window. 



A. Dynamical analysis : adiabatic approach 

Using a complete adiabatic approach for our system 
implies the assumption that the large spin's movement 
is slower compared to the electrons jumping through the 
system. This is an additional assumption on top of the 
mean- field approach, where quantum fluctuations are ne- 
glected already. But even for the derivation of the rate 
equations one needs an adiabatic approximation for the 
electrons tunneling into the system. We performed the 
latter in the last section to decouple the lead-transition 
functions, cf. Eq. (fT3|) . A full adiabatic approach also 
assumes that the electronic spin changes on a time-scale 
which is much smaller than the one of the large spin. This 
approximation for the interaction between the electronic 
spin and the large spin reduces the number of dynam- 
ical equations to three, because solely the ones for the 
large spin remain. The equations of motion for the elec- 
tronic spin operators are solved with the help of Green's 
functions. 

We can use the adiabatic approach to search for fixed 
points and also perform a rough characterization of them. 
The predictions for the classification of the fixed points 
coincide quite good with the actually obtained nonadia- 
batic results. However, a complete adiabatic treatment 
of the system cannot catch the whole dynamics of the 
system. Even in the infinite bias regime the adiabatic 
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results do not coincide with the ones obtained from the 
rate equation or the nonadiabatic approach. The rea- 
son therefore lays in the lack of additional friction terms, 
which could be gained by expanding the Green's func- 
tions as done in RefJ^ for a related system. 



1. Fixed point analysis of V : center or saddle point 

The fixed points , where the large spin is completely 
polarized parallel to the magnetic field, cf. Eq. ([T8]) . ap- 
pears also for the adiabatic system. But the correspond- 
ing value of the z-component for the electronic spin does 
now depend on several system parameters; for zero tem- 
perature it yields 



<5;(±j)> = Efe arctan 

a 

Tc4 



2irT 



arctan 



Oc* - ft) 

r/2 



r/2 



(22) 



The x- and ^/-component are zero as for the rate equa- 
tion approach. For infinite bias the ^-component solely 
depends on the tunneling rates (S* lB (±j)) = (IYt — 
r u )/2r, which coincides with Eq. (ffgj) . 

For the fixed points , one eigenvalue of the Jacobian 
is zero and the two other read 



£ 0l± - ± 7 



\ 



( 



d(S x ) 



V 



d(J x ) 

with T = X(S* z (±j)) + B. 



\ 



(23) 



Therewith two possible realizations can appear. If the 
root is real the eigenvalues are complex conjugate and 
the fixed point is therewith stable and can be classified 
as a stable center. This means, that periodic oscillations 
for the adiabatic spin dynamics are expected. The real 
part is equal to zero and therefore no damped oscillations 
should appear. This missing damping is characteristic 
for the adiabatic results, hence this approach cannot de- 
scribe the complete dynamics of the system. This is valid 
even in the infinite bias regime, since the appearing dy- 
namical features like self-sustained oscillations require an 
additional positive and negative damping^. 

In Fig.[7J (S*) as a function of the applied bias is de- 
picted, together with the real (gray dashed-dotted lines) 
and the imaginary (cyan solid line) part of the corre- 
sponding eigenvalue. The two upper graphs show the 
behavior in regime I for (J*) polarized parallel (left) and 
anti-parallel (right) to the magnetic field. We observe 
regions where the imaginary part of the eigenvalues is 
equal to zero and therewith the oscillations disappear. 

The ranges of finite real part are slightly different for a 
parallel and a anti-parallel polarization of the large spin 
z-component. These differences appear for small bias, 
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FIG. 7: Results of the dynamical analysis for region I and 
II. Here, parameters are chosen as in Sec. lII Bl The solid 
black line denotes the ^-component of the electronic spin 
(S z ) and the external bias is chosen in a symmetric man- 
ner /xl = —Mr = Vbias/2. Zero imaginary part appears 
in region I for parallel alignment of (J z ) in the range of 
Vbias/2A G [—5; 60] and for anti-parallel alignment the range 
yield Vbi as /2A £ [10; 60]. For region II the imaginary part 
disappears in the range of Vbias/2 A £ [— 5; oo] for the parallel 
state and Vbi as /2A £ [6.4; oo] for the anti-parallel state. 



the intervals of finite real part are denoted in the caption 
of Fig. [TJ For Vq the intervals with finite real part are 
smaller as for Vq . Depending of the existence of other 
fixed points, the regions where only one fixed point has 
eigenvalues with finite imaginary parts, are promising for 
the occurrence of spin-flips. 

Based on the chosen polarization of the leads and 
therewith the different tunneling rates for the spin-down 
electrons we obtain a system which is not symmetric. 
Hence, the dynamical behavior for negative detuning is 
different than for positive detuning V^ias > 0- This is 
clearly visible in Fig.[7J where for large negative bias the 
imaginary part in region I is always unequal to zero. The 
infinite bias result for region I equals (S*) IB = ±0.25 
as expected, whereby the sign depends on the detuning. 
For region I the system performs periodic oscillations for 
a bias Vbias/2 A > 60. 

In contrast, in region II (2nd row in Fig. [7)) no finite 
imaginary part for the high bias regime exists. This co- 
incides with the infinite bias result for the nonadiabatic 
approach, where the concerned fixed points do not ap- 
pear as stable centers or spirals. But here, for a small 
bias range and also for negative detuning, we find, that 
the eigenvalues can become complex and therewith the 
fixed points can exist as center. For positive detuning 
the fixed point 7^ has already turned into a saddle point, 
but the anti-parallel state is alive until Vbi as /2A = 6.4. 

For region III, which means for B z /\ > 0.25, the be- 
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FIG. 8: Results for the fixed points Vg N in region I: numer- 
ically calculated density plots for (S z ) + B z /X and (S x ) as a 
function of the large spin components. These functions must 
be both equal to zero for the appearance of V± N . The yellow 
(blue) area in the first two columns corresponds to positive 
(negative) values. Along the borders of these areas the func- 
tions are equal to zero. The right column depicts the super- 
position of the first two columns, there a contact of the blue 
and yellow area denotes a solution for V^ N (red circles) . The 
bias increases from the upper to the lower row, explicit values 
are Vbi as /2A = 5; 8; 10. The numerical results were sorted by 
their sign, with the purpose of highlighting the positive and 
the negative areas. The bifurcation appears approximately at 
V bias /2A & 7. 



havior of the fixed points is similar to region I, if the 
magnetic field and the tunneling rate are small. But for 
higher values of these parameters the real part of the 
eigenvalues is always zero. 



2. Two conditions for the fixed points V 



± 

SN 



The next fixed points obtained from Eq. (j6j) include 
the requirement that (S x ) = and the condition (S z ) + 
B z /X = has to be complied with as well. We define 



SN 



: ({j*),±^f-(j^-{j*)Mr z )) 



(24) 



In the infinite bias regime V^ N coincide with the four 
fixed points in region II Eq. (pjjj), obtained from the rate 
equation. The investigation of V$ N in the finite bias 
regime is possible only numerically. In Fig. [8] we present 
results obtained for region I. The cycle shapes originate 
from the conservation of the large spin, hence the radius 
equals j. Outside of these cycles (J y ) becomes complex 
and no real solution for V^ N exists. 

The stability is defined by the eigenvalues of V^ N . One 
eigenvalue of V^ N is again zero, due to the fact, that 
we actually deal with a two-dimensional system. The 




FIG. 9: Results of the dynamical analysis of Vg X for region 
I and II. Outside of the depicted range, the fixed point is 
complex. The real range is limited by (J y ), which is only real 
inside of the dashed line. Parameter as in Sec. lII Bl 



remaining eigenvalues yield 



c 2,3 — 2^ y ' 



d(S x ) d(S z ) 
d(J z ) d(J x ) 



± 



d(S x ) , d(S z ) 



Ms x )d(s z ) 



\ \d(J z ) d(J x )J d(J x )d(J z ) 



(25) 



For a small bias, Vbi as /2A < 5, only one solution for 
Vg N exists, corresponding to the case (J*) = 0, which 
we name as Vg X , the latter is stable until T4ias/2A < 7 
cf. Fig.O By further increasing the bias another bifur- 
cation appears, where V sx becomes unstable and two 
other fixed points are created. These points are symmet- 
ric to the (J x ) - axis and move with higher bias values 
further to (J z ) = ±j. When they reach the border of the 
cycle, they disappear and (S z ) + Bz/X = has no longer 
a real solution. The emerging fixed points have nega- 
tive/positive real eigenvalues and can be characterized 
as stable/unstable nodes. For region II these points do 
not disappear for high bias and the nodes are the only 
remaining stable solutions of the system. These nodes 
correspond to fast damping behavior for the spin compo- 
nents. 

For V^ x , where (J x ) = 0, the last term in Eq. ([25]) is 
unequal to zero and the component (J*) is obtained from 
the transcendental equation 



{I^larctan 




\27rr 


r/2 



Tc4 

27TT 



arctan 



r/2 



(26) 



with £| ^ = Ed ± 0.5(B Z + A(J*)). This equation has no 
solution for (J*) in the infinite bias case and therefore 
the fixed points Vg X ^° not ex i s ^ there. In the regime 



of a low magnetic field, the right side of Eq. (|26|) has 
to be small. Following from that, the argument of the 
arctan function has to be small as well, leading to the 
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FIG. 10: Fixed points for unpolarized leads as a function 
of tunneling rate Y and magnetic field B z . The symmetric 
bias is chosen as Vbi as /2A = 5. The right graphs depict the 
adiabatic results for the large spin, the small cycles on the 
top/bottom correspond to and the larger cycle to 
(region A). Magnetic field for region A/B is B z /X = 0.1/0.3. 



estimate for the evolution of the z-component as linear 
to the applied bias (J*) ~ T^ias- 

This linear behavior is clearly visible in Fig.[9j where 
we plotted V^x anc ^ ^ s eigenvalues as a function of the 
applied bias. The arcs corresponds to the (J*) compo- 
nents, which determines whether the fixed point exists 
or not. Within this arc the point is real and therewith 
physically reasonable. The radius of the arc is limited by 
the length of the large spin. 

For both regions we observe small ranges with finite 
imaginary part, where the fixed point can be classified as 
a stable center. As well as ranges, where a saddle point 
occurs. In regime I, the point T$x starts its existence 
for Vt,ias/2A ~ — 5. At the bias value, Vq turns into a 
saddle point, cf. Fig. [71 

This behavior matches quite good those of Vq. We 



can interpret Vg N as the complementary points in the 
region, where disappears. As we discussed before, 
for a certain bias region the stable solution turns into 
a saddle point and we can estimate, that then the fixed 
points appearing in Fig.[5J correspond to stable solutions 
of the dynamical system. 



3. Fixed points for unpolarized leads 

For unpolarized leads the tunneling rates become spin- 
independent T = Tl + Tr. There only the fixed points 
Vq and Pg X exist, as depicted in Fig.lTOl In dependence 
of the tunneling rate and the magnetic field, two main 
regions appear. The first one, region A, has the stable 
fixed points Vq and Pf x , these points can be character- 
ized as centers. In this region spin-flips of the large spin 
can be possible, because Vq is no stable solution. But 
the system is still multistable and therewith the trajec- 
tories can also end up in the fixed points P$ x . Note, 



that for an total adiabatic ansatz, a complete spin- flip is 
not observable, because the trajectories starting parallel 
to the magnetic field end up in P^ x , as shown in the 
right graphs of Fig.lTDl If one reverses the direction of 
the magnetic field, region A contains Vq instead of Vq . 

In region B solely Vq are stable fixed points. By fur- 
ther increasing the bias, the region where Pf x exist gets 
smaller and in the infinite bias case only the second region 
persists and all electronic spin components (Si) become 
zero. Than the two systems decouple and the large spin 
oscillates with the Larmor frequency. 



B. Results for the nonadiabatic approach in the 
finite bias regime 

We expect to obtain the same fixed points for the nona- 
diabatic approach as for the adiabatic approach. For the 
case Vq , we want to show, that this fixed point also 
appears in the nonadiabatic regime. There, all lead- 
transition functions for different spin B acr > vanish and 
we obtain for the ones with equal spins the simple result 



= i-^G?M, (27) 

Z7T (uj - E% + iky) 27T 



containing the spin-dependent single-level Green's func- 
tion Gp(uj) without coupling of the two electronic levels, 
due to (J*) = 0. Also the equation for the z-component 
of the electronic spin is straightforwardly obtained from 



a a. 



arctan 



(Ma - e%) 



r/2 



(28) 



This result coincides with Eq. ([22]) and the fixed point 
is identical to the one obtained with the adiabatic ap- 
proach. In the same manner, we can construct also the 
other fixed points. In the general stationary case, all 
lead-transition functions decouple from the equation sys- 
tem and can be expressed by retarded Green's function. 
The adiabatic expressions for electronic spin components 
can be reconstructed, too. 

This equivalence is not valid considering the adia- 
batic eigenvalues, they are not directly transferable to 
our nonadiabatic approach. The fixed points are long- 
time quantities and stationary solutions of the dynamical 
system, in contrast the eigenvalues include higher order 
terms and their calculation requires the first derivations 
of the system's variables. 

But some predictions derived from the adiabatic eigen- 
values also appear in the nonadiabatic approach. There- 
fore, we can estimate that the behavior for the full time- 
dependent solution is strongly influenced by the adiabatic 
eigenvalues. 
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FIG. 11: Nonadiabatic results obtained in region I for three 
different bias values, which is chosen in a symmetric manner 
VWs = Ml/2 = —fiji/2. The electronic current, depicted 
in the lowest row, is separated into left/right and spin-up/- 
down contributions. The magnetic field yield B z /X = 0.1 and 
r/A = 9. Initial conditions as in Sec. lII Bl for regime I. 
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FIG. 12: Nonadiabatic results for region II. The upper row 
corresponds to a small bias Vbi as /2A = 2, where the spin com- 
ponents oscillate around the coordinates of Vq . Increasing 
the bias, leads to a disappearing of the oscillations and the 
fixed point V^ N is present, see lower row for Vt>ias/2A = 5. 
The magnetic field yield B z /X = 0.2 and r/A = 0.7. Note, 
that the time-scale of the current and the spin operators differ 
slightly. Initial conditions as in Sec. lII Bl for regime I. 



1. Region I: disappearing of the oscillations 

We start our presentation and discussion of the nona- 
diabatic results with the ones obtained in region I, corre- 
sponding to a low value of the external magnetic field and 
a large tunneling rate. In the case of infinite bias, we ob- 
tained an oscillating behavior between two cycles. This 
is different to the rate equation result, where damped 
oscillations appear. 

From the adiabatic analysis of the last section, we ob- 
tained the prediction, that the fixed point Vq of the sys- 
tem turns from a center into a saddle point for the bias 
range Vbias/2A G [—5; 60]. When Vq loses its stability at 
^bias/2A w — 5 a supercritical pitchfork bifurcation^ — 
appears and V^ x are created, they exist until Vbi as /2A « 
7 as centers. If the latter fixed points become unstable, 
two stable nodes are born and we expect the oscillations 
to disappear. Note, that these predictions originate from 
the adiabatic eigenvalues. 

As depicted in Fig.QTJ we obtain these features also 
in our nonadiabatic results. There, the spin compo- 
nents and the current are depicted for three different 
bias values. For ^i as /2A = 2 we observe non-sinusoidal 
self-sustained oscillations. (S z ) oscillates around —B z /\, 
which corresponds to V^ x . 

If we increase the applied bias for the initial conditions 
corresponding to the results depicted in Fig.QTj the os- 
cillations disappear in the range of ^ias/2A ~ 8. They 
run into one fixed point of the kind Vf N , which is clearly 
visible in Fig. El were (S x ) = and (S z ) = -B z /\. 

By varying the initial conditions, we observe that not 
all oscillations have disappeared for a bias in the range 
of H»ias/2A « 8. Choosing the initial conditions near to 
the anti-parallel fixed point Vq , we observe oscillations, 



which are smoothly sinusoidal and run around Vq . For 
^bias/2A w 10, the oscillations disappear and the tra- 
jectories run in the same fixed point as depicted in the 
middle graph of Fig.QTJ This result coincides perfectly 
with the border predicted in the adiabatic analysis. 

The bifurcation point, where a revival of the oscilla- 
tions appear, is also correctly predicted by the adiabatic 
analysis. For Vbias/2A « 60 we observe again oscillations 
around Vq . These oscillations slightly show the two cycle 
behavior as in the infinite bias case, but the differences 
between the radii of the cycles is not large. The latter is 
visible in the third column of Fig.[TTJ where the results 
for Vbias/2A = 100 are depicted. 

The lower row in Fig.[TTJ depicts the electronic current 
(X aa (t)), split of into its constituent parts. The right 
tunneling amplitude is equal to zero and following from 
that, as well the corresponding current channel. There- 
fore, the current 2r^ should be equal to the total current 
through the system, hence current conservation remains 
valid. For the nonadiabatic current, this is ensured for 
the time- aver aged current values, but can be different in 
the time-dependent case. There, we observed some accu- 
mulation of current in the central region. 

When we consider the current channels for Fbias/2A = 
2 in Fig.[TTJ we notice that all channels reach their 
minima if the (J z ) component is maximal and in this 
case close to j. For this low bias regime, the shift of 
the energy level, due to the coupling to the large spin 
±1/2(B Z \ + (Jz)), leads to a positioning of them slightly 
outside of the transport window. Therewith the current 
is minimal. 

The current corresponding to the spin-up electrons 
flows in the direction of the bias. For the spin-down elec- 
trons the situation is more complicated, because they are 
not allowed to leave the system through the right lead 




0.5 



FIG. 13: The large spin's ^-component for several bias values 
in regime II. Para meter as for Fig.ll2land initial conditions are 
(X) = (J y ) = V50-9.9V2, (J z ) = 9.9 and (S y ) = (S z ) = 
0, (S x ) = 0.5. 



due to = 0. The electrons can either flip their spin, 
stay in the lower level or flow back into the left lead. If 
the last process appears the electron moves against the 
bias and the current becomes negative. This feature is 
slightly visible in the lowest graph of the first column in 
Fig. [HI there the current Z14 drops below zero for a quite 
small region. 

Again, we can address this effect to the position of the 
large spin's ^-component. In the case Z14 < 0, the shift- 
ing of the energy level is the other way around, because 
(J z ) is negative. Following from that, the spin-down level 
lays above the spin-down level and additionally in the 
neighborhood of the left Fermi edge (/xl = 2/A) and elec- 
trons can occupies empty states in the left lead. After 
(J z ) passed its minima, the spin-down level moves down 
and therewith its current channel drops as well. 



2. Region II: Negative Detuning and spin-flip of the large 

spin 

The feature of negative current is better visible in re- 
gion II, as we see in the current results for ^ias/2A = 2, 
which are depicted in the upper row of Fig.[l2j There, 
the large spin's z-component oscillates close to —j and 
therewith, both effective levels are clearly situated out- 
side the transport window and the spin-down level lays 
again above the spin- up level. The oscillations of (J z ) 
are comparatively small and do not influence the behav- 
ior of the effective levels as much as for region I. They 
stay outside of the transport window for all times. 

We try to interpret the evolution in time for the current 
channels focusing on transition between the levels. If Zjx 
is negative, the right current reaches its minimum and 
the left current for spin-up electrons is maximal. Due to 
< ZfL? we can assume that spin- flips from the lower 
to the upper level happen and a depletion of the upper 
level into the left lead appears, due to the negativity of 
2|l- This is in accordance with the evolution for the 
(S z ), which decreases in this range point. 

The other way around, when (S z ) increases, the left 
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FIG. 14: The graphs depict the different contributions to the 
current in regime III. The right graph presents a detail of 
the left graphs drawn together for a short time range. Initial 
conditions and parameters as in Sec. Ill Bl for regime III. 



current for spin-down electrons is maximal, as well as 
2^r. But ZfL is close to zero in this regime. Therefore, 
we interpret this kind of current cycle in the following 
way. A spin-down electron enters the upper level, flips 
into the lower lead, due to the interaction with the large 
spin, and finally leaves the central system to the right 
lead. This would explain why the current for spin-up 
electrons is maximal at the right lead and minimal at 
the left lead. 

If the bias is increased, the regions of negative spin- 
down current vanish. As well as the oscillations of the 
spin components, as depicted in the lower row of Fig.[T2l 
Remember, for the infinite bias regime, we found no peri- 
odic oscillations at all in this region. In accordance with 
the foregoing adiabatic analysis, the trajectories enter 
one fixed point of kind V SN . 

The adiabatic analysis, also predicts an earlier disap- 
pearing of the fixed point Vq , corresponding to a paral- 
lel alignment of the large spin, as the anti-parallel fixed 
point Vq . From that we propose the possibility of spin- 
flips for the large spin in the region, where only Vq has 
finite imaginary part: ^ias/2A G [—5; 6.4]. 

To monitor how the first fixed point disappears we 
choose our initial conditions close to Vq and assume 
negative detuning. The results are depicted in Fig. [131 
For Vbias/2A = — 10 the trajectories perform smooth os- 
cillations around Vq . As expected from Fig.[7J the z- 
component of the electronic spin (S z ) ~ 0.2 and is there- 
with positive. As discussed in the last section, the ob- 
served system is not symmetric and therewith the fixed 
point Vq stays alive by further decreasing the bias. 

Increasing the bias we see a different behavior, 
(J z ) drops down when the bias passes the threshold 
^bias/2 A = —5 as it is clearly visible in Fig. [131 At first 
the trajectory runs into the fixed point Vg X , laying in the 
middle of both points . But even for 14ias/2A = —4, 
the spin components enter the anti-parallel fixed point 

If we decrease the tunneling rate, we observe a tran- 
sition to chaotic-like oscillations as in the infinite bias 
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FIG. 15: Stable fixed point as a function of tunneling rate and magnetic field. The external bias is increased from left to right, 
the explicit values are denoted in the graphs (in units of [A]). Additionally, the infinite bias result is plotted in the last graph 
on the right side. The blue area denotes the regime where no oscillations appear and only the nodes from the 'P 5iV -series exist. 



case, see Fig.|4j For small values of the bias, the chaotic- 
like behavior is a little suppressed and the trajectories 
oscillate comparatively smooth with a high frequency. 



3. Region III: oscillations and high frequency 

In Sec. Ill Bl we presented results for larger values of 
the magnetic field, where the trajectories oscillate around 
. This kind of behavior is recovered for the finite bias 
regime. But we observe slight differences in the regime of 
small bias. In Fig.[]3]the current for regime III is depicted 
(M)ias/2A = 2), the results for the left and the right cur- 
rent are split into their contributions from spin-up and 
down electrons. We observe that the current channel Z14 
oscillates around the zero axis and following from that, 
the current is flowing in both direction. The frequency 
of the oscillations is quite high and the spin components 
oscillate between two cycles, which signatures are clearly 
visible in the left graphs in Fig.[T4l 



IV. CONCLUSION 

With our nonadiabatic approach we are able to probe 
the rate equation and the adiabatic approach and thus 
can discuss the advantages and disadvantages of these 
two methods. 

The rate equation approach fromi^, is a quite practical 
method. The numerical effort is low and the dynamical 
system can be investigated analytically. However, this 
method neglects higher order transitions and we learnt 
from the nonadiabatic results, that these are relevant in 
the regime of low magnetic field. In the latter region, 
the rate equation misses parts of the nonlinear dynam- 
ics. If the magnetic field increases, the results of both 
approaches coincide. 

The limitation of the rate equation approach to high 
external bias is another disadvantage: at finite bias, we 



observe even richer dynamics within the nonadiabatic ap- 
proach. There, features like spin-flips of the large spin 
and suppression of the oscillations as a function of the 
applied bias appear. The disadvantages belonging to the 
nonadiabatic approach are the high numerical effort and 
its inaccessibility for further analytic investigations. 

In our work, we utilized the adiabatic Green's func- 
tion method to analyze our system. However, we omit- 
ted the presentation of time-dependent results, because 
the complete adiabatic approach does not capture a lot 
of the dynamics even in the infinite bias regime. Chaotic, 
damped or two cycle-like oscillations are not observable 
within this approach: the dynamical system reduces to 
three equations of motion, which is not suitable to de- 
scribe the system in the range of low magnetic fields. 
Only if the system performs smoothed low-frequency os- 
cillations or runs directly in a stable fixed point, the adi- 
abatic approach can correctly describe the dynamics. 

On the other hand, the adiabatic Green's functions are 
suitable for the analysis of the system. For the finite bias 
regime, our initial adiabatic analysis matched quite well 
the behavior of the nonadiabatic approach. Therewith 
we can re-define the dynamical regions as depicted in 
Fig. [151 It is clearly visible that the external bias is an 
important parameter for this system. The regions where 
only one of the main fixed points exists are promis- 
ing for switching of the large spin, but due to the fact 
that we deal with a multistable system, there are always 
additional stationary solutions present. 
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Appendix A: Frequency doubling of (S z (t)) 

The equation of motion for the z-component of the 
electronic spin in Eq. JT3J) does not directly couple to the 
magnetic field, but it couples to the product of the x- 
component for the large spin and the ^-component of 
the electronic spin. The latter develop sinusoidal in time 
with the same frequency uo s and without a phase shift, 
cf. Fig. [5] If we approximate their evolution in time with 
(J%ff(t)) = ajs'mu s t and (Syff(t)) = a e s'muj s t, we can 
estimate an effective solution for the z-component of the 
electronic spin (A s = aja e ) 



(Sf"(t->oo)> 



Lt 



■ 14 



■2A S ) 



2r 

cos(2cj s t) + 2u s /T sin(2cj s t) 



El) 

4 I 



(Al) 



Hence, the oscillation goes with twice the frequency of 
the other spin components. Note, that for a more general 
ansatz, p.e. (J%?^ if)) — a j sincj s t + bj cosc<j s t, the result 
is similar and differs solely in the prefactor. 



Appendix B: Adiabatic approach: Green's functions 

The expectation values of the spin operators in fre- 
quency space are obtained from 



(Si) (a;) = -- tr [G < ( W )(7 i ] , i = x,y,z, (Bl) 



containing the Pauli spin matrices 



1 

1 



-i 

1 



1 
-1 



If we reconsider the effective Hamiltonian Eq. ([7j) for this 
system, we remember the similarity to a parallel two-level 



system and as in that case, the Green's functions have 
matrix character 



G(w) 



G n (w) G t ;(w)\ 
G#(u>) G u (u>)) ■ 



The off-diagonal functions refer to the coupling to the 
x-component of the large spin. Without this coupling we 
end up with two independent levels. For the derivation 
of the electronic spin we need the lesser Green's function 
whose calculation takes usage of the Keldysh equation 

G< a , = ]T G^M <#,,„,( W ), (B2) 



with the lesser self energy 

a a 

The retarded and the advanced Green's function are ob- 
tained from their equation of motion and yield 



ou — e a — z^o-fj 



1(4) 



UJ ■ 



fW 2 

(B4) 



with the retarded/ advanced self energy 



E 



\V k0 



OJ - Skacr ± iO 



( B5 ) 



which real part A aa (uj) solely leads to a shift of the level 
energies Sd and therefore can be neglected. The imag- 
inary part corresponds to the tunneling rate T aa (uj) = 
r 

-L (y rr • 
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